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Abstract 

This paper offers analytical capability for 
aerodynamic parametric studies and sensitivity analysis of 
rotary wings in axial flight by using a 3 -dimensional un- 
distorted wake model in curved lifting-line theory. The 
governing equations are solved by both the Multhopp In- 
terpolation technique and the Vortex Lattice method. The 
singularity from the bound vortices is eliminated through 
the Hadamard’s finite part concept Good numerical 
agreement between both analytical methods and finite dif- 
ference methods are found. Parametric studies were made 
to assess the effects of several shape variables on 
aerodynamic loads. It is found, for example, that a rotor 
blade with out-of-plane and inplane curvature can 
theoretically increase lift in the inboard and outboard 
regions respectively without introducing an additional in- 
duced drag. 


Nutation 

AR aspect ratio, R / c t 

a lift curve slope, =27t 

b number of blades 

c chord length 

c 0 chord length at blade root 

C| chord length at blade tip 

c e equivalent chord 

c c(x)lc 0 

C L total lift coefficient 

C D total induced drag coefficient 

i 

C p total induced power coefficient 

1 r 

K coefficient factor, c t / c a 

M number of terms in Multhopp Interpolation 

Technique or number of panels in VLM 
m summation integer, also weighting factor 

index 

p tip shape parameter 

R rotor radius 
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radial distance 

area of a rotor in X-Y plane, =R c t 
normalized coordinates with respect to R 
tip length parameter 
twist shape parameter; also tarf '(eT|) 
inplane curvature 
twist angle 

0W/e o 

twist angle at blade root 
twist angle at blade tip 
out-of-plane curvature 

dummy variable representing spanwise 
distance 

taper ratio, c x !c 0 
air density 

rotor solidity, bcJnR 

twist ratio, 9|/8 0 

inflow ratio, = V, + 0^/9^ 

climb inflow ratio, V cluni /R Q 

induced inflow ratio, Uf+Oj- 

induced inflow ratio from curved lifting-line 

induced inflow ratio from trailing vortex 

lines 

tan - '(ex) 

azimuthal angle; also Digamma function 
circulation ; also Gamma function 

2 

normalized circulation, TlR Q0 O 
angular velocity 

Introduction 


A key ingredient in developing reliable and efficient 
procedures for design optimization of aerospace structures 
is sensitivity analysis. This is especially true for 
aerodynamic sensitivity analysis which is of growing im- 
portance for aerodynamic shape optimization as well as 
the interdisciplinary process of integrating aerodynamics 
with other disciplines, A recent paper by Sobieski stres- 
ses the need for aerodynamic sensitivity analysis and 



makes a plea to the aerodynamic community for extend- 
ing their present capability to include sensitivity analysis. 
The present paper attempts to partially fill this void for 
the rotary wings. 

In recent years there has been a resurgence of interest 
in the idea of reducing the induced drag of an aircraft by 
wing design. Some examples of recent work suggest the 
use of inplane 2 * 4 or out-of-plane 5 curvature. There have 
been few attempts to extend this idea to rotary wings. The 
present paper is intended to establish sensitivity analysis 
for rotor blades in axial flight. 

In previous work, Rand and Rosen 6 developed a lift- 
ing line model which is capable of calculating the 
aerodynamic load along inplanc curved-blades of a 
helicopter in hovering and axial flight Their method, 
however, should be classified as the Vortex Lattice 
Method (VLM) and is an approximate approach to solve 
the singular integro-differcntial equation. The present 
paper deals with the same problem by using an analytical 
approach to solve the more general governing equation 
which can include the analysis of out-of-plane curved 
blades. The analytical approach, in the present paper, uses 
the Multhopp Interpolation technique which is a special 
case of Fourier Series method and widely used in fixed- 
wing cases for both lifting-line and lifting-surface ’ 
models. As mentioned in Ref. 6, the difficulty, raised in 
the case of inplane curved blades, is the fact that the 
velocity which is induced on the lifting line by the bound 
vortex itself obtains infinite values. This phenomenon is 
due to the singularity of the order 1/Ixl which occurs in 
the development. This singularity may be handled 
elegantly by the Hadamard’s finite part concept 10 as 
described in Ref. 11. In this paper, we apply a similar 
technique to eliminate the singularity. 

This paper also contains parametric studies which 
can quickly assess the effects of several shape variables 
on aerodynamic loads. The design variables are taper 
ratio, tip length, tip shape, aspect ratio, inplanc curvature, 
out-of-plane curvature, twist ratio, and twist shape. Most 
rotor blades may be described by these shape parameters. 
Aerodynamic functions arc local circulation, local in- 
duced velocity, total lift coefficient, total induced drag 
coefficient, lift-drag ratio, and induced power coefficient. 
The rotor performance is strongly dependent on these 
aerodynamic functions. 

Aerodynamic Functions and Shape Design Variables. 

The present paper deals with a rotor system having b 
identical, equally spaced blades. These blades rotate about 
the shaft with a constant angular velocity Cl.. The rotor in 
Figs. 1 and 2, has a radius R and no hinge offset The 
steady lifting-line theory, based on the three-dimensional 
undislortcd wake model, will be applied to cases for 


which helicopter rotors or propellers are in axial flight 
(including hovering). The shape design variables charac- 
terize three essential parameters of a blade: lifting-line 
shape, chord distribution, and twist distribution. The 
design variables are taper ratio, tip length, tip shape, 
aspect ratio, inplanc curvature, out-of-plane curvature, 
twist ratio, and twist shape. Aerodynamic functions are 
local circulation, local induced velocity, total lift coeffi- 
cient, total induced drag coefficient, lift-drag ratio, and 
induced power coefficient Further details of the shape 
design variables are given next. 

Lifting line shape 

The reference blade in Figs. 1 and 2 is positioned at 
the 90 degree azimuthal angle. The origin or center of 
rotation is on the blade root one quarter-chord away from 
the leading edge. The inplane position y, and out-of-plane 
position z, of the lifting-line are assumed to have the 
form: 

7=y(jc) = « 2 ; Z = g(x) = KX 2 

where x is the spanwise direction, £ is the inplane curva- 
ture, and k is the oul-of-planc curvature. All quantities 
arc normalized with respect to blade radius, R, When e, 
and k are equal to zero, it is a straight lifting-line. The 
rotor blade can have either backward inplane curvature 
(e > 0) or forward inplane curvature (e < 0) and have 
either upward curvature (k<0) or downward curvature 
(k > 0), Also, we assume that lei, and lid are smaller than 
0.3 so that the lifting-line theory remains valid. 

Chord distribution: 

The chord length c(x) normalized with respect to root 
chord has the following variation along the spanwise 
direction : 

c(x) = &L = [ 1 + x(\ - 1 )] 1 1 - jr ,/p ] p 
C o 

where X is the taper ratio and is defined as the ratio of tip 
chord c x to root chord c py shown in Fig. 2. Note that we 

define X as the inverse of conventional taper ratio. When 
X is equal to 0, the blade is triangular. 

The tip shape parameter is denoted p and is always 
positive. The physical meaning of p is to define the tip 
shape of a blade as shown in Fig. 3. The tip length 
parameter is denoted P and is also positive. The physical 
meaning of P is to define the extent of the tip taper as 
shown in Fig. 4. In general, all cases of chord distribution 
can be defined on the combinations of these three 
parameters. 

Another important design variable is the solidity a, 
(or the aspect ratio A/?,) which is defined as: 

_K = _b_ 

nR ~n(AR) 


2 


where b is the number of blades, AR is the aspect ratio, 
and c e is the equivalent chord, which is based on the 

equivalent thrust generated by the reference rotor with a 
rectangular planform, and has the following form: 

/ xTc(x)dx 

c t = = K(mX$p)c 0 

\xTdx 

0 

where m is the weighting factor index and is dependent 
on design criteria, such as blade area or area, thrust, or 
torque weighted solidity. Once m, (J, p and X are specified 
for a rotor, K can be analytically expressed in terms of the 
Gamma function as : 

K = K x + (\-l )K 2 
where 

_ (m+l)pr(p+l)r(p^ ffl P) 

1 n/H-l+p + mP) 

(m+i)proH-i)r(2P + /»B) 

2 T(>f 1+2P + mp) 

In order to calculate sensitivity derivatives, it is 
necessary to derive 9 K/d X, 9 K/d p, and 9 K / 9 p, which 
axe analytically expressed in terms of the Gamma and 
Digamma functions as r and y respectively: 



= X'JyO^ 1) - \|/(/h- 1+P + mP)l 
+ (X- i)/f 2 ry(p+i)-y(>+i+2p + mP)i 

()k r i 

% = Kl P + ( m+ 1 ) “ ( OT+ 1 ) ¥(/”■ 1 +P + xnP) 

+ (X- 1)^| + (m+ 2) y(2p+mp) - (m+2) y(p+ l+2p + mP) 

Twist distribution: 

The normalized twist angle of attack 0 (jc) with 
respect to twist angle at blade root along the spanwisc 
direction is defined in the following form: 

= i + , 8 <t-i) 

where x is the twist ratio, defined as the ratio of tip twist 
0! to root twist 0 O . 6 is the twist shape parameter and is 

always positive. The physical meaning of x is to define 
the twist distribution along a blade as shown in Fig. 5. 
The reason behind this x definition is that an ideally 
twisted blade can be approximately modelled by concave 
twist (i.e., (kxcl). In order to calculate sensitivity 
derivatives, it is necessary to derive 3 0/3 x, and 3 0/3 5 
which arc expressed as: 
dO(x) _ s 

9t x 


= ln(jc)/(T-l) 

Theoretical Derivations 


Governing equation for circulation: 

To derive the governing equation for circulation, we 
assume that the flow field is incompressible. Each blade 
is represented by a curved lifting line. The circulation 
along each lifting line is T(r). Since F varies along the 
blade, vortex filaments must trail behind the rotating 
blade by dT/dr . It is further assumed that the trailing 
vortex system is an undistorted wake and moves 
downward with a climb velocity, V climb , which is equal to 

the sum of the rotor axial velocity and the averaged in- 
duced-velocity around the rotor disk. The derivation fol- 
lows similar derivations which have been developed in 
Refs. 12-13. However, an additional term of downwash 
due to the self-induced velocity from the bound vortex of 
lifting line has to be included. A similar study for fixed- 
wing cases can be found in Ref. 14. The^ governing equa- 
tion used to solve for local circulation, T, is very similar 
to those contained in other references. See for example, 
Ref. 6, Eq. (34) or Ref. 15, Eq. (3.17). It is the integro- 
differential equation which has the following form: 

^c(x)[x0(x) - D(x)l = ^T(x) 
where 


w(jc) = t),W = v L (x) + Ujix); 


be, Kbc. _ 

r(jt): 






nAR nR tiR ' 
here, a is the lift slope =2 n and v c is the climb inflow 
ratios. x> r and are the induced inflow ratios, con- 
tributed by the curved lifting-line and trailing vortex lines, 
and can be derived by Biot-Savart Laws as: 




U*) 


= — j 1 


4n °bc — t|I[(e+k) 2 (jc + t|) 2 + 1] 

dT 




/ \ ”*1 f 1 v dl 

”' <X) = 4*1 

r tl 2 (l + 6 2 TlVx(l + E 2 X 2 ) Q5 Tl(l + £ Z Tl 2 ) 0 - 5 CO S (vi/ + 5-0 

0 |/?t 3 

where 

- 2x(l+e 2 x 2 ) a5 Ti( 1 +e 2 Ti 2 ) n5 cos(y+S-0+(w c y-KTi 2 +)cx 2 ) 2 
and 

y = y + 2 /” C = tan _1 (£x); 5 = tan -1 (er|) 
b 


3 



Note that is the induced velocity due to the bound vor- 
tex of an inplane curved blade. We neglect the effect of 
the curved lifting-line from other blades. However, the in- 
fluence of the trailing vortex lines from other blades has 
been included. 


Equations_for aerodynamic function: 

Once r(jc-) and o/jc-) are obtained along the spanwise 

position jc-, (by solution of governing equation) the total 
lift coefficient C L total induced drag coefficient Cn> and 


total induced power coefficient C P are: 


& 


pr£lT(r)dr 


ARjsnxpAx 

*0 p (ROfSfl. J 

C D f^rQw{r)T{r)dr 


0„ 


p (RQ)%B 0 
C, 


= AR YjXj r (x^-otyAXj 


C P C, C D 

TT ■ ^ 17 


j j 

where S R is the blade area and is equal to Rc e . (Note that 
our definition of S R is different from those normally used 

in the helicopter field.) If we divide the equation of total 
induced power coefficient C P by C D , it gives an index to 

judge how good a rotor is by the use of lift-drag ratio, 
C L /C D . This design index is slightly different from those 

in the fixed-wing field 2 " 5 , which use ct/c D - 


Governing equations for derivatives of circulation: 

To obtain aerodynamic derivatives, we use the con- 
cept developed by Yates in Ref, 16. Take partial deriva- 
tive of Eq. (1) with respect to a typical design variable 
denoted Q as: 


3 Jga- 


c(jc)[x0(x) - X)(x)] 


} - 


( 2 ) 

BQ\2K l 

Solve Eq. (2) for ihc derivative of circulation, 31730. The 
governing equations for the derivative of circulation with 


K AR 3F 3v 
kc 3X + 3X 
KAR 3f dv 
nc 3p 3P 

KAR 3f t dv 
lie 3 p 3 p 


KARj 


Jtcfl+^l+X)] 


3 i 

(K AR 

3X1 

1 nc 


rcep 2 


3P 


KAR 


K AR lnO-y^ 1 ) 3 

dp 


nc p 2 


TtC 

KAR 


nc 


f (3) 
- (4) 


- (5) 


KAR 3r ( 3\) -K r 

(6) 

Ttc 3 AR dAR nc 


KAR 3r + 3u .30 

(7) 

nc 3x 3x 3 x 


KAR dr 3\) 30 

(8) 

nc 35 + 35 *35 


KAR dr 3\> 
tic 3k + 3k 

(9) 

K AR dr dx> „ 
+ — = 0 

(10) 

7tc 3e 3e 



where 3 v>/3 Q, in Eqs. (3)-(8), is expressed as: 

b _ 

3t>(x) _ -1 rf( 3 T/3 Q) p 
3 Q ~ 4jr*o i ^ dx\ -*o 

r| 2 ( 1 +z 2 r\ 2 )-xT\^( 1+e V)( l+£ 2 Ti 2 rcos(y^S-0 ( 1 

am 

However, 3 o/3 k and 3 *0/3 e have the following forms: 
b _ 

3 o(jc) _ —I f * rf( 3 T/3 k) 

3k 4^0" ch\ " 


/i-l 



[Y— 

AK^o^dr] 

n= 1 


(q 2 -xqcos(H/+S-0)(q 2 -x 2 )(\) c v-Kii 2 +tcr 2 ) 


c 

b _ 

3o(.x) _ -1 r*y d(3 T/3 e) 


dydr\ 


( 12 ) 


d n 


r q 2 ( 1 +e 2 q 2 )-xq Vp+£ 2 x 2 )( 1 +eV) cos(\j?+5-0 . , 

o l*t 3 




£j ^(Aq 2 +2eVq 2 )cos(v+^0+xn(e^-lj(x-q)sin(v+S-0 ]i . 

V(l+e 2 jc 2 )(l+e 2 q 2 ) 

3 f 1 y 3f r r\ 2 ( 1 +e 2 Ti 2 hxq V( l+eV^ 1 +e 2 ~n 2 ) +cos(y+5-p 

dnJ n 


4n J o^di] J o 


r 4 4 

\ZT[ +ZX - 

mi (x 2 +Ti 2 +2£V)cos^+&-0+^n(£ 2 ^- 1 x*— 0 )sin(\j7+5-Q ^ ^ 
V(l+e 2 x 2 )(l+mn 


4 


t X-j' 1 * < x + 1 f ( e + 1 ~ 3 M E ' l ~1> f^) d T, <13) 

n °_Lc - • ql [( e + K) 2 (jr + T |) 2 + 1] ' 

Once d r (x t )/d Q and d u ( (jc,)/3 Q are obtained along 
the spanwise position or,, (by solution of governing equa- 
tion) and the d C L /d Q, d C D /d Q, and dC p /d Q can be 

easily found through the chain rules from equations of 
C L C D , and C P . 

Numerical Solutions 

Solutions for circulation: 

Two numerical methods have been developed to 
solve Eq. (1) for circulation. These are (a) the Multhopp 
Interpolation technique, which is a special case of Fourier 
Series method and widely used in fixed-wing cases and 
(b) the Vortex Lattice Method, which can be easily ex- 
tended for lifting-surface theory. 

Multhopp Interpolation Technique : The idea is to 
transform the spanwise coordinate by : 
r = l-cos(q)) 

2 

and replaces the downwash integral \)^(jc) in Eq. (1) with 
midpoint trapezoidal rule summation to obtain: 




-fi- 

4k j q , dr\ 


n=\ 


r 

¥ A 


ti 2 ( 1 +eV)-;cTiV(l +eV)( 1 +e 2 q 2 ) cos(\j/+5-C,) 


l^l 3 


M , v 


= — X- 

4M^|d<p 


dych\ 
(14) 


where 


- iC 


n=l 


q 2 ( 1 +e 2 q 2 )-jc ( Ti,V ( i +e 2 r 2 )( 1 -t-e 2 !) 2 ) cos(\g+8 r Q ( 1 4 ' a ) 

rfy 


' 2 = ti 2 (l+e 2 q 2 )+x 2 (l+e 2 r 2 ) 

- 2jc,ti V( 1 +e 2 r 2 )( 1 +c 2 q 2 ) cos(y+5-Q+(t> c y-KT| 2 +Kjt 2 ) 2 
and 

|^| 2 = il^l+e 2 !!?^- ( 1 +e 2 x 2 ) 

- 2jc ( t| / V( 1 +e 2 xf)( cos(\j/+6 r C i H(\) c v-)cn 2 +K^) 2 

Equation (14) representsUie downwash due to M trailing 
vortices of strength (<f T/J (p )^( n/M), which are along 

the helical wake and are located at: 


R- 


(/ - 0.5) 7t , , 

<P, = ’ ,= I -- W 

with control points located at: 

The Multhopp Interpolation formula suggests that: 

- 2 V-V (15) 

r (<P) = Z r jL sin(m<p>) sin(mcp) 

j= 1 m=l 

where T - represents T((p / ) and 

Equation (15) is based on the following orthogonality 
property; 

M - 1 

2 

— £ sinfmepj) sin(w(p ( ) = 5 y 

m=l 

here, 8^ is the Kronecker-delta . When Eq. (15) is used 
in Eq. (14), o 7 >(x ( ) reduces to : 

A/ M - 1 A/-1 

u 7< x «) * 773 Z Z Olm sin(ffi<p:) cos(m(p ; ) 

2Af 


A/-1 


/=! y=l iw=l 


F J 

j=i 


(16) 


where 


M M-X 


aJj = "j XX™ sin(/mp-) cosfm^) P l7 

2M 1= 1 m=l 
T 

here, - are the influence coefficients of trailing vortices 
and can be solved for by numerical methods without any 
difficulty, if x i and ti / are specified and P tl can be calcu- 
lated. However, P u contains an integral with an infinite 

upper limit and no closed-form solution exists. A numeri- 
cal solution has been developed by the current authors in 
Ref. 17 to calculate P tl scmi-analytically. Results are very 

accurate with minimal computation time. 

Similarly, if we replace the downwash integral \) L (. x ) 
in Eq. (1) with the Multhopp Interpolation formula at 
control points x i% then x> L (x-) becomes: 

W = H — r — r? f ^ dr \ 

n °Lc 4 — T|l[(e+K) 2 (jc t - + q) 2 + 1] 


T^aI 


1 nr\)di\ 


4re J ol jc 4 - — T) 1 

Af-l M - 1 it 

£ Try sin(m<p ) sin ( m y) sin (9) <L <P 

Z, l jZ, sin l m( P/> J | cos((p ) _ C0S ( (p ) | 


2Mn 


=1 


w-i 


H 


(17) 
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where 



M~\ n 

vH 1.****P J 

m= I 0 


sin(m(p) sin((p) d <p 
I cos(<p ( ) - cos(<p) I 


here, A~ are the influence coefficients of bound vortices 

and contains an integral with a Log singularity. This in- 
tegral diverges when <p- is equal to <p. It should be inter- 


preted in other than the usual sense. The answer is given 
by the concept of the "finite part" established by 
Hadamard 10 and has been applied to fixed-wing cases in 

Ref. 11. In Appendix A of this paper, we show that A- 
can be obtained analytically if T is expressed as the 
Fourier sine series. In Eq. (17), we apply the assumption 

of small (k -f e) 2 . Cases for large (k + e) 2 can also be cal- 
culated elegantly and are described in Appendix A. 

Finally, substitution of Eqs. (16) and (17) for the 
downwash in Eq. (1) and application of resulting equation 
at spanwise location <p- gives a (A/“l)x(Af~l) matrix 


equation for Af-1 values of T's as: 

tf-i _ _ (18) 

X A i] r j - x i 0 « - v 0 o . * = 1 ,~M-l 

j = 1 
where 

A^A^Ajj + ^KAR/inCi) 

Once V's are obtained, C L , C D ,C L I C D , and C P , can 

j is. is. r - 

be calculated as described before. The local induced 
velocity t>- (jc-) is found as: 

Af-l 

= V*,) + ^d x d * X ( A fj + A Jj ) 0 

Vortex Lattice Method (VLM): The conventional 
VLM discretizes the rotor span into a number of panels, 
M. The trailing point is located at the boundary of each 
panel, and the control point is placed at the midpoint of 
each panel However, there is a controversy as how to 
choose collocation points. Several papers * * ’ have 
shown that this conventional VLM converges slowly and 
docs not appear to approach the correct limit for fixed- 
wing cases. A similar phenomenon has been found for 
rotary-wing cases in Ref. 15. AUfficsc papers suggest that 
the semi-circle collocation points, which are defined in 
the Multhopp Interpolation technique, arc essentially re- 
quired for VLM when the fastest rate of convergence is 
desired. Most researchers have selected collocation points 
arbitrarily or equally-spaccdly. Such selection may be 
adequate for aerodynamic analysis. But, it is not suitable 
for aerodynamic sensitivity analysis. 

The VLM in this paper uses the semi circle colloca- 
tion points to solve Eq. (1). For the trailing downwash 
integral i) 7 (x) in Eq. (1), we apply the first order finite 


difference on T with midpoint trapezoidal rule summation 

T 

to obtain A- as: 


4 = i (J’u-'V i) 


(19) 


« 4n v 0 

where P i} is defined by the Eq. (14-a). A common way to 
calculate A ^ is described in Ref. 6. The straight vortex 


filament with unit circulation is stretched from each trail- 
ing point. Then, the sum of this vortex filament and Eq. 
(19) is equivalent to the influence of a horseshoe vortex 
at panel j at control polnTx-. However, this approach is 

not completely rigorous, since the curved lifting-line has 
been approximated by several straight vortex filaments. 
Another way (used in this paper) to solve the problem is 
to apply Eqs. (17) and (19) for downwash calculation in 
Eq. (18), since both equations have the same control 
points. 


Solutions for derivatives of circulation: 

Two numerical approaches can be used to solve Eq. 
(2) for derivatives of circulation and result in same 
numerical solutions. The first one is to discretize Eq. (2) 
by following the path of discretizations used for solutions 
of circulation. The second approach simply differentiates 
a set of equations (expressed in equation number (18)) 
with respect to a particular design variable Q , Then, we 
solve a matrix equation for A/-1 values of derivatives of 
circulation, which can be expressed as: 

( 20 ) 

a e rj.i-iju-i 


v\ d _h v'£4 

X A ij a q dQ X 

H -H 

where 


dA it a4 3 aJj a f K AR/( jc c ( ) ] 

dQ = dQ + dQ + 8, > dQ 
Equation (20) is a set of linear equations and can be 
solved if the right hand side of Eq. (20) is known. How- 
ever, not all terms exist for all cases. For example, the 
term 3 A*j/d Q exists for nonzero inplane curvature e; the 
term 3 ( x t 8- )/3 Q is valid for parameters of twist dis- 
tribution; and the term 3 [ K AR/{ nc i )]/dQ exists only 


for parameters of chord distribution. The procedures are 
systematically analogous to analyses for circulation which 
are described in the previous section. 


Resulls 


Test problems 

Equations (1) and (2) have been normalized with 
respect to some reference quantities. The following 
parametric values will be used in this paper for the refer- 
ence-blade: 8 0 =1, X=l, p = |5 = 0.5, \> c = 0.06, 

5^0.6, K = z = x-m=O i AK=10 {or a = 0.0637); 
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b- 2. In some cases, we compare results with dilTcrcni 
chord distributions which arc defined as: 

(1) elliptic blades: X = I,p = 0.5, p = 0.5 

(2) rectangular blades: A. = 1, p = 0, P = 0 

(3) triangular blades: X = 0,p = 0, P = 0 

(4) tapered blades: X = VS, p - 0, P = 0 

(5) parabolic blades: X = 1 , p = 0.5, P = 0.5 

where all blades have the same solidity based on the same 
blade area (i.e., m = 0). Figure 6 shows the normalized 
quarter chord distribution along the spanwise direction for 
these blades. 

Effect of number of panels or terms, M 

Table 1 is the convergence results of circulation (T) 
and derivative (3173k) with respect to out-of-plane curva- 
ture (k) for an elliptic blade at midspan x = 0.5. This table 
implies that aerodynamic functions and derivatives are 
equally accurate from both Mullhopp Interpolation for- 
mula and VLM but are sensitive to the number of dis- 
cretization panels and terms M. Also, the analytical 
methods are confirmed by the finite difference method 
which solves Eq. (1) twice with k = 0 and k = 0.001. The 
similar accuracy is also obtained for C L C D ,C P and its 

derivatives. For a reasonable number of terms or panels 
(M=20), the methods provide approximately 4 digits of 
accuracy for aerodynamic functions and 2 digits of ac- 
curacy for aerodynamic derivatives. The following results 
for this paper arc based on M=20 by the Multhopp Inter- 
polation formula. 

Effect of curvatures of the the lifting line 

Figures 7-9 are results of C L ,C D and C L /C D as 

I i 

functions as out-of-plane curvature k for various blade 
shapes. Figure 7 shows that C L is nearly proportional to 

k. However, C D is nearly constant for each blade shape 

i 

shown in Fig. 8. Results of C L /C 0 in Fig. 9 have the 

similar trend as those for C L . Therefore, by curving the 
blades in the opposite direction of thrust (i. e., k > 0), the 
rotor would be more efficient The reason is because the 
trailing vortex lines move away from the blades when k is 
increased. Consequently, the induced velocity along the 
span is decreased and the thrust is increased. This 
phenomenon is clearly shown in Figs. 10 and 11. Figure 
10 presents the variation of circulation along the span. 
Figure 1 1 presents the variation of induced velocity. It is 
interesting to note that the inboard region (x = 0 -» 0.8) is 
improved when k is increased. Unfortunately, such a 
blade is physically unrealizable due to the fact that the 
helicopter blades have to cone up (i. e., k<0) to compen- 
sate the centrifugal forces for structural considerations. 


For effects of inplane curvature e, similar trends of 
results have been found as shown in Fig. 12 for CJC D . 

The rotor lift-drag ratio will improve when e is increased. 
Therefore, by sweeping the blades backward (i. e., e > 0), 
the rotor would be more efficient Figure 13 presents the 
variation of circulation along the span. Figure 14 presents 
the variation of induced velocity. Note from Figs. 13 and 
14 that the outboard region (x = 0.4 -» 1.0) has been im- 
proved when e is increased. 

One might wonder how much additional improve- 
ment might be made by the combinations of k and e. One 
answer is shown in Fig. 15 where contour lines of CJC D 

i 

arc plotted for different values of k and e. It clearly 
shows that positive k and e can always improve the lift- 
drag ratio, C L /C D , of a straight blade (i.e., k = e = 0). For 

l 

example, numerical results indicate that making 
k = e = 0.3 can increase CJC D by 16 percent as com- 

i 

pared to the straight blade. However, making k = e = -0.3 
can decrease Lhis ratio by 10.4 percent. Also, Fig. 15 indi- 
cates that the effect of k is slightly stronger than the ef- 
fect of £. Figure 16 presents the contour lines of 
3 {CJC D )f d £ which represent the slope of Fig. 15 with 

respect to £. The peak value occurs around £ = 0 and 
k = 0.05. It is observed, from Figs. 10 and 13, that dual 
peak circulations appear. This is very common for rotors 
which do not have an ideal twist. Unlike the conventional 
blade element theory, the three dimensional undistorted 
wake model in this paper can catch this "dual peak" 
phenomenon. Cases for different twist distributions and 
inflow ratios have also been investigated. All results indi- 
cate that positive k and e are very promising parameters 
to increase lift in the inboard and outboard regions 
respectively without introducing an additional induced 
drag. 

Effect of chord distribution 

Figures 7,8,9, and 12 indicate that elliptic and 
parabolic blades are better than rectangular blades. For 
example, parabolic blades always provide more lift with 
the same induced drag than rectangular blades. Elliptic 
blades can provide the maximum C L and C D . For all 

I 

blade shapes, a case with larger C ^ is always accom- 
panied with larger c D . and vice versa. Although triangular 

i 

blades provide the maximum C L /C D , they are not 

i 

desirable. The reason is because C L and C D are also very 

i 

small and this design is unrealistic due to the need for 
larger C L to provide sufficient lift. A compromise design 

is by the use of tapered blades, as shown in Figs. 7 and 9. 
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Such blades can provide the minimum desired lift with 
reasonable lift-drag ratio. 

The results for effect of aspect ratio are shown in Fig. 
17-20. Figure 17 indicates that C L grows with increased 

aspect ratios for all blades. However, C D , shown in Fig. 

18, increases when aspect ratio is smaller than 9. Then, it 
decreases with a larger aspect ratio where the flow field 
is approaching the two-dimensional case and the induced 
velocity tends to decrease along the span. Figure 19 
shows the results of C L IC D and indicates that C L /C D is 

nearly proportional to aspect ratio. Figure 20 is the result 
for the derivative of C L /C D: Increasing the aspect ratio is 

increasing the lift-drag ratio. However, the resulting in- 
crease in blade weight would tend to limit the increase in 
aspect ratio. Results for variations of p y P , and X and its 
derivatives have also obtained and will be presented for 
CJC D only. Figure 21 indicates that variation of ( -d C D i 

with tip length parameter, p, for an elliptic blade has the 

minimum value around P ~ 0.25. This minimum value is 

around p ~ 0.20 with tip shape parameter, p, as shown in 

Fig. 22. Figure 23 shows that variation of C L lC D with 

* 

taper ratio, X , for an elliptic blade has the maximum value 
at X = 0.0 which represents the triangular shape. 

Effect of twist distribution 

The results for effect of x are shown in Figs. 24-26. 
Figure 24 indicates that C L is nearly proportional to x. 

This is due to the fact that a large value of x provides 
large angles of twist along the span. It is interesting to 
note that, from Fig. 25, decreases when x is negative. 

Consequently, An optimum design of a rotor can be 
specified around -0.05 <x< 0.05, which are shown in 
Fig. 26 with maximum values of C L jC D . The results for 

the effect of 8 are shown in Fig. 27 which presents the 

variation of C L iC D . It seems to imply that an optimum 

* 

design is around 0.35 < 8 < 0.40. However, This is only 
partially true due to the fact that a small value of x always 
provides small values of C L and 

Conclusions 

In summary, this paper offers analytical capability for 
aerodynamic parametric studies and sensitivity analysis of 
rotary wings in axial flight by using a 3-dimcnsional un- 
distorted wake model in curved lifting-line theory. 

Analytical expressions are developed for aerodynamic 
functions and sensitivity derivatives with respect to shape 
design variables including taper ratio, tip length, tip 
shape, aspect ratio, inplanc curvature, out-of-plane curva- 


ture, twist ratio, and twist shape. The governing equations 
of aerodynamic functions and derivatives are solved by 
using both Multhopp Interpolation Technique and Vortex 
Lattice Method. Numerical results show that aerodynamic 
functions and derivatives are insensitive to discretization 
methods but are sensitive to the number of discretization 
panels. For example, for a reasonable number of terms or 
panels (M=20), the methods provide 4 digits of accuracy 
for aerodynamic functions and 2 digits of accuracy for 
aerodynamic derivatives. 

Also, it is found that a rotor blade with out-of-plane 
and inplane curvature can theoretically increase lift in the 
inboard and outboard regions respectively without intro- 
ducing an additional induced drag. However, positive out- 
of-planc curvature is physically undesirable due to struc- 
tural considerations, while positive inplane curvature z is 
easier to obtain. 


A ppendix A 


In Eq. (17), the finite part ( FP ) for the following in- 
tegral which contains a 1/lxl type singularity is defined in 
Ref. 11, Appen. A, as follows: 


FP 


J ^ , d q 

J 0l*i - T|l 


lim 


0 U 0 
di\ + 2^x-)lne] 


d H 


(A-l) 


+ f' 

If the variable in Eq. (A-l) is transformed from r\ to <p 
according to the semi-circle points, An equivalent equa- 
tion can be obtained as follows: 
n 


/ = FP j 


M sin((p) d (p 
I cos(q> ( ) - cos(<p) I 


lim [J 


c ^-^0 


7((p ) sin(<p) dj p 
I cosfrp,) - cos(tp) I 


f 1 M sin (<P ) d <P 
I cos(<p ( ) - cos(<p) I 


+ 2/(<p,) Ine ] 


(A-2) 


where 


f® ■ 

e = sirup 

Equation (A-2) can be solved analytically if /Tn) can 
be expressed as the Fourier sine series. For example, let 
/(ti) be equal to sin Then, 

/ = fpl sin (wcf>) sin(<p ) d (p 

* I nr\c( \ I 


I cos(<p ( ) - cos(cp) I 


1 


= 2 [ • /w_1 ) - C (<P , ^+i) ] 


(A-3) 


where 


GOP, jn) = FP J «£ i gfl 4 L 

‘ J I cos(«p ( ) - cos(tp) I 
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The integral G is obtained from the recurrence relation: 

„ 4 sin (mq> ) 

G(<p- ,m+ 1) - 2 cos ( <p-) G(<p- ,/n) + G(<p- ,m-l) = 

tn 


where the initial integrals are: 

G(<p ' >0) = ii^) ,n(jl>,((p ‘ )) 

G(<p, ,1 ) = 2 (p— re + 4 cot (sin ( <p ( )) In (sin ( <p ( >) 
Substituting Eq. (A-3) in Eq. (17) wc obtain 

= ~~ X sin(/n<p ) j fin(m<p> sin<q>) cp 
" ; J I cos(<p ( ) - cos(<p) I 


2Mk 
M - 1 


m=l 


0 


4A/tc 


X sin(m<fy) [ G(<p, zn-1)- G(q>, jn+\) 


(A-4) 




Equation (A-4) is only good for small (K + e) 2 . For cases 

for IargcJ»c + e) , we use Multhopp interpolation formula 
for /i(q) F(ti) as : 

Af-l A/-1 
2 

A r = — £( AHy X sin(m(p y ) sin(mtp) 


where 
A(q) = 


h 1 


1 

1 


[(£+K) 2 (x. + T1) 2 + l] 1 * 5 


1 


(A-5) 


1.5 


[(e+K) 2 (1 - 0.5cos(cp ; ) - 0.5cos((p)) 2 + 1] 
Substituting Eq. (A-5) in Eq. (17) we obtain 

= 5 

4 M n [(e+K ) 2 (1 - 0.5cos(<p ( ) - 0.5cos(q> ; )) 2 + 1] * 


Af-l 


X sin(m<p y ) [ G((p- jn- 1) - G(<p, ^n+1) ] 


m=I 


This work was performed at the NASA Langley Re- 
search Center under NASA contract NAS 1-19000 for 
Lockheed Engineering and Science Company, technical 
monitor Dr. Howard M. Adclman. 

References 

'Sobicski, J., "The case for Aerodynamic Sensitivity 
Analysis," Sensitivity Analysis in Engineering, NASA 
CP-2457, 1986, pp. 77-96. 

2 Hcndcrson, W. P., and Holmes B. J., "Induced Drag 
- Historical Perspective," SAE Acrotect ’89 Conference, 
Anaheim, California, September 25-28, 1989. 

3 Burkett, C. W., "Reductions in Induced Drag by the 
use of Aft Swept Wing Tips," Aeronautical Journal . 
December 1989, pp. 400405. 


4 Gucrmond J., " A Generalized Lifting-Line Theory 
for curved and Swept Wings," J. Fluid Mech .. 1990, vol. 
211,^.497-513. 

Lowson, M., L., "Minimum Induced Drag for Wings 
with Spanwise Camber," Journal of Aircraft. Vol. 27, NO. 
7, July, 1990, pp. 627-631. 

Kand, O. and Rosen A., "A Lifting Line Theory for 
Curved Helicopter Blades in Hovering and Axial Flight," 
Jour nal of AHS. vol. 29, no. 3, July 1984, pp. 3-10. 

'Dcjameltc, F. R„ "Arrangement of Vortex Lattice on 
Subsonic Wings," Vortex-Lattice Utilization, NASA SP- 
405, 1976, pp. 301-323. 

8 Multhopp, H., "Methods for Calculating the Lift 
Distribution of Wings," British A.R.C., R&M 2884, 
January, 1950. 

9 Lamar, J. E., "A Modified Multhopp Approach for 
Predicting Lifting Pressures and Camber Shape for Com- 
posite Planforms in Subsonic Flow," NASA TN D-4427, 
July, 1968. 

°Hadamard, J., Lectures on Cauchy’s Problem in 
Linear Differential Equation . Dover, 1932, pp. 133-153. 

1 'Gibbs, C. P„ Lifting Surface Theory for Swept and 
Crank Subsonic Wings . M. S. Thesis, North Carolina 
Stale University, Raleigh, 1985. 

12 Baskin, V. E„ et al„ Theory of the Lifting 
AilSfpi, NASA TTF-823, 1976, pp. 147-163. 

Rand, O. and Rosen A., "The Aerodynamics of 
Helicopter Rotors in Hovering and Axial Flight Using 
Lifting-Line Models," TAE No. 485, TECHNION, Haifa, 
Israel 1982. 

Ashenbcrg, J. and Weihs, D., "Curved Lifting-Line 
Theory for Thin Planar Wines." Israel Journal of Technol- 
Ogy, Vol. 20, 1982, pp. 160-165. 

duccd-Flow Calculations bv Optimum Choice of Mesh . 
Ph. D. Thesis, Georgia Institute of Technology, Atlanta, 
1988. 

Yates, E. C., Jr., "Aerodynamic Sensitivity from 
Subsonic, Sonic, and Supersonic Unsteady, Nonplanar 
Lifting-Surface Theory," NASA TM 100502, 1987. 

1 Chiu, Y. D. and Peters, D. A., "Numerical Solu- 
tions of Induced Velocities by Semi-Infinite Tip Vortex 
Lines," Journal of Aircraft . Vol. 25, NO. 8, August, 1988, 

pp. 684-694. 

18 

Lan, C„ "A Quasi-Vortex-Lattice Method in Thin 
Wing Theory,” Journal of Aircraft . Vol. 11, No. 9, 1974, 
pp. 518-527. 

19 

Guermond J., "Collocation Methods and Lifting- 
Surfaces," Eur. J. Mech.. B /Fluids . 8, No. 4, 1989, pp. 
283-305. 


9 


Table 1 Convergence of circulation and derivatives 
with respect to k for an elliptic blade at midspan. 


Number 
of terms/ 
panels, M 

circulation at 
midspan, f (0.5) 

circulation derivatives at midspan, 3 f/3 k 

Analytical Methods 

Finite Difference Method 

Multhopp 

VLM 

Multhopp 

VLM 

Multhopp 

VLM 

8 

0.01950 

0.01970 

0.00468 

0.00470 

0.00468 

0.00470 

16 

0.01951 

0.01957 

0.00468 

0.00468 

0.00467 

0.00469 

24 

0.01951 

0.01953 

0.00470 

0.00470 

0.00470 

0.00470 

32 

0.01951 

0.01952 

0.00471 

0.00471 

0.00471 

0.00471 

40 

0.01951 

0.01951 

0.00472 

0.00472 

0.00472 

0.00472 

48 

0.01951 

0.01951 

0.00473 

0.00473 

0.00473 

0.00473 

56 

0.01951 

0.01951 

0.00474 

0.00474 

0.00474 

0.00474 
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Fig. 2 Top view of rotor blade. 


lJ r quarter chord 



x 


Fig. 3 Physical meaning of dp shape parameter, p, with 

(3 = 0.01 ; X = 1. 


Fig. 5 Physical meaning of twist shape parameter, 5, with 

t = 0. 



x 


Fig. 6 Test cases for blades. 
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Fig. 7 Variation of C L with out-of-plane curvature, k , for Fig. 10 Radial distributions of r(x) of elliptic blades for 
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Fig. 8 Variation of with out-of-plane curvature, k, 
for various blade shapes. 


Fig. 1 1 Radial distributions of i) ( {x) of elliptic blades for 


various k. 
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Fig. 9 Variation of C L /C D with out-of-plane curvature, 
ic, for various blade shapes. 


Fig. 12 Variation of C L /C D with inplane curvature, e, 
for various blade shapes. 
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Fig. 13 Radial distributions of r(r) of elliptic blades for 


various e. 


Fig. 16 Contour lines of derivatives of C L /C D with 
respect to t with k and e for an elliptic blade. 
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Fig. 14 Radial distributions of o t {jc) of elliptic blades for Fig. 17 Variation of C L with aspect ratio, A*, for various 

various €. blade shapes. 



22o r c 0 


190 " /// 


* II. TIT* 

J TS.W31 
7 ntw 
i itnu 


->* ELLIPTIC 
-B- 7APER<1/3) 


RECTANGULAR TRIANGULAR 


Fig. 15 Contour lines of C L /C D with k and e for an 
elliptic blade. 


Fig. 18 Variation of C D with aspect ratio, AR , for 
various blade shapes. 
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Fig. 19 Variation of C L /C D with aspect ratio, AR, for 
various blade shapes. 



4 8 8 ^ 10 12 14 is 

Fig. 20 Variation of derivative of C L /C D with aspect 

I 

ratio, AR, for various blade shapes. 
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Fig. 21 Variation of C l /C D with tip length parameter, 
P, for an elliptic blade. 


Fig. 22 Variation of C L /C D with tip shape parameter, p , 
for an elliptic blade. 
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Fig. 23 Variation of with taper ratio, X , for an elliptic 
blade. 



-0.3 -0.2 -0.1 0 0.1 0.2 0.3 

X 

Fig. 24 Variation of C L with twist ratio, x, for various 
blade shapes. 






Fig. 25 Variation of Cp with twist ratio, t, for various 
blade shapes. 


Fig. 27 Variation of C L /C D with twist shape parameter, 

i 

8, for various blade shapes. 



Fig. 26 Variation of C L /C D with twist ratio, T, for 
various blade shapes. 
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